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In this work, we investigate calcium titanate (CaTiOs - CTO) using X-ray diffraction and Raman 
spectroscopy up to 60 and 55 GPa respectively. Both experiments show that the orthorhombic 
Pnma structure remains stable up to the highest pressures measured, in contradiction to ab-initio 
predictions. A fit of the compression data with a second-order Birch-Murnaghan equation of state 
yields a bulk modulus Kq of 181.0(6) GPa. The orthorhombic distortion is found to increase slightly 
with pressure, in agreement with previous experiments at lower pressures and the general rules for 
the evolution of perovskites under pressure. High-pressure polarized Raman spectra also enable 
us to clarify the Raman mode assignment of CTO and identify the modes corresponding to rigid 
rotation of the octahedra, A-cation shifts and Ti-O bond stretching. The Raman signature is then 
discussed in terms of compression mechanisms. 
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I. INTRODUCTION 



In the past, investigations of ABO3 perovskite-type 
oxides have been a rich source for the understanding of 
structural mechanisms not only in perovskites but also 
in oxides and other ionic materials in general. Beyond 
their fundamental interest, perovskites attract continu- 
ing attention as functional oxides, often related to dielec- 
tric, ferroelectric and multiferroic properties [J, 0|- The 
ideal cubic structure of ABO3 perovskites is rather sim- 
ple, with corner-linked BOg octahedra and the A cations 
sitting in the space between the octahedra 0-111 • Most 
perovskites present, however, structural distortions away 
from this parent cubic structure, which can be driven by 
external parameters such as temperature, pressure, stress 
etc. One of the major challenges in perovskite science is 
to identify these distortions, which are often rather subtle 
and thus difficult to follow. 

Here we will focus on a class of materials in which the 
structural distortion is dominated by tilts (rotations) of 
the BOq octahedra Q. This kind of distortion is ferroe- 
lastic and corresponds to anti-ferrodistortive (AFD) in- 
stabilities at the zone boundary It is well known [H 
that such a tilt distortion (the tilt angle) can be driven 
to a large extent by temperature or pressure as exem- 
plified in the model materials SrTiO^ 0, il or LaAlOg 
[9l-[ll| . It is generally accepted that temperature reduces 
tilt instabilities, i.e. the tilt angle decreases with increas- 
ing temperature. The situation is less straightforward for 



the effect of high-pressure: based on a pioneering work 
by Samara et al. [l^ . it was for a long time believed 
that pressure increases systematically AFD tilt instabili- 
ties, i.e. the tilt angle increases with increasing pressure. 
However, more recently, this view has been challenged by 
a number of counter-examples (lol. ITsl. IT^ . This situation 
has motivated the search for criteria able to predict this 
evolution with pressure. Generally speaking, the volume 
reduction in perovskites can be accomodated by two dif- 
ferent mechanisms [isl [l5l [l6| : bond compression that 
favours a reduction of the distortion, and tilting of the 
octahedra that results in an increase of the distortion. 
The pressure dependence of the distortion is then con- 
trolled by the relative compressibilities (3 of the BOq and 
AO12 polyhedra (or equivalently of the A-0 and B-0 
bonds): Pb/Pa < 1 implies an increase of the tilt angles 
under pressure and vice versa. Experimentally, Andrault 
and Poirier [l7| first estimated the individual polyhedron 
compressibilities using the empirical relation between the 
volume and the bulk modulus of the polyhedra given by 
Hazen and Finger [l^. More recently Zhao, Ross and 
Angel [15I . [l6| have developed a model based on the cal- 
culation of the valence bond sums, whose main feature 
is the prediction that the tilt angle should increase un- 
der pressure for A^+B^+Os perovskites but decrease for 
B^~^03 compounds, such as aluminates AAIO3, in 
agreement with experimental data. Moreover, the bond 
valence sum model provides a way to estimate the angle 
change under pressure [l6| . 

Having the evidence that octahedra tilts can either in- 



2 



(a) 



(c) 




3 1 




020 


"^■^1 23 
1,22 

' 222 






'^^1' 212 


f 








■ 311 

10 1 




i 


2 2 

1-11 


0-2 




= 2-12 






■ 3-1 

2-2 1 



Figure 1: (a) Polarized light micrograph of a slice cut from 
the as-grown crystal, showing large ferroelastic domains, (b) 
Picture of the pressure chamber at 60 GPa during the X-ray 
diffraction experiment, (c) Example of a diffraction pattern 
integrated over the full angular range with peak indexing. 



crease or decrease with increasing pressure, it is now in- 
teresting to ask if there are perovskites in which the tilt 
angle does not change with pressure, but where the com- 
pression is largely dominated by bond compression, and 
up to what point such structures can remain stable. To 
the best of our knowledge, the search for such materi- 
als has until now received only little attention, with the 
notable exceptions of CaTiOg (CTO) and CaGeOg [Hi 
whose tilt angles and spontaneous strains show only lit- 
tle variations up to 10 GPa. CTO has been in the past 
often investigated both at low temperatures for its in- 
triguing dielectric properties [1^ and at high pressures 
and temperatures for its suggested analogy with geologi- 
cally relevant perovskites such as MgSiOa [2l| . Previous 
investigations at high pressure have been carried out us- 
ing Raman spectroscopy up to 26 GPa j 2 1 | and X-ray 
diffraction in various pressure ranges |l9l |22|. [23| . Al- 
though Xiong et al. |22| have claimed a transition to a 
hexagonal phase at approximately 10 GPa, this has not 
been confirmed by either the Raman scattering experi- 
ment by Gillet et al. [2l|, or the diffraction experiment 
by Wu et al. 

In this work we carry out experiments on CTO sin- 
gle crystals at room temperature (RT) by using Raman 
scattering and X-ray diffraction. The general aim of our 
study is to follow and interpret the structural distortions 
up to high pressure, and more specificall y: ( i) to clarify 
the presence of a theoretically predicted [2J| and exper- 
imentally disputed (2l] - l23| pressure-induced phase tran- 
sition and (ii) to explore the compression mechanisms at 
very high pressure up to 60 GPa. 



II. EXPERIMENTAL DETAILS 

CaTiOa single crystals were grown by the floating-zone 
technique, with powders of CaCOs and Ti02 (99.9% pu- 
rity) as starting materials. The as-grown crystals showed 



a clear domain structure visible under polarized light, 
with approximately 100 |xm sized domains (figure [IJa)). 
The samples used for the experiments were polished to 
a thickness of about 10 |xm with a lateral extension of 
10 to 30 iJ.m and did not show any domain structure, 
neither under polarized light nor in the diffraction pat- 
terns. All experiments were performed in diamond-anvil 
cells (DAC) with various pressure transmitting media, as 
detailed in the following. The cells have the Boehler- 
Almax diamond design with a cuUet of 300 |J.m. The 
pressure chamber was sealed by a stainless steel gasket 
pre-indented to a thickness of about 40 |J.m. We checked 
by measuring the gasket thickness after the experiment 
that the pressure chamber was always larger than the 
crystal in order to rule out any bridging of the diamonds 
by the sample. A picture of the pressure chamber at 60 
GPa is shown in figure [IJb). 

X-ray diffraction experiments were performed on the 
ID27 beamline at the ESRF. The beam was monochro- 
matic with a wavelength of 0.3738 A selected by an iodine 
K-edge filter and focused to a beam size of about 3 |^m. 
The signal was collected in the rotating crystal geometry 
(-30 ° < a; < 30 ° ) on a marCCD detector. The crystal 
has its c axis close to the DAC axis, as can be infered from 
the indexed pattern given in figure [TJc) . A precise cal- 
ibration of the detector parameters was performed with 
a reference silicon powder. Helium was used as a pres- 
sure transmitting medium. The fluorescence of ruby was 
used as a pressure gauge j25| . During our experiments, 
we have carefully followed the splitting between the two 
fluorescence lines Rl and R2 of ruby: the splitting was 
approximately 1.40 nm at ambient conditions and drifted 
slightly under pressure up to a maximum of 1.50 nm at 
the highest pressure measured. 

Raman scattering experiments were performed with 
neon as a pressure transmitting medium. The spectra 
were recorded on a Labram spectrometer with a low fre- 
quency cutoff at 100 cm^^. The exciting laser line was 
the Argon at 514.5 nm. The laser power was kept at 10 
niW on the DAC to avoid heating of the sample. 



III. X-RAY DIFFRACTION 

At ambient conditions, CTO is orthorhombic with 
space group Pnma [2^ . The oxygen octahedra are tilted 
along the three directions of the cubic perovskite unit cell, 
with a tilt system a~6+a^ using Glazer's notation 
The presence of twins is rather common in Pnma per- 
ovskites in general and may complicate considerably the 
analysis of the diffraction pattern. In our case however, 
the analysis of the superstructure reflections correspond- 
ing to the tilts revealed that our crystal was a single 
domain crystal, with only very weak reflections arising 
from a possible twin. The diffraction spots could there- 
fore be indexed unambiguously in the Pnma space group. 
The diffraction pattern did not show any major change 
over the entire pressure range investigated and the same 
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Table I: Lattice parameters as a function of pressure, given as 
pseudo-cubic lattices parameters related to the orthorhombic 
parameters tto, feo and Co by ape = ao/\/2, bpc = feo/2 and 
c„c = C0/V2. The pressure can be considered hydrostatic up 



to 40 GPa only (see text). 



indexation was used throughout the experiment. The 
lattice parameters were determined by a least squares fit 
to the positions of the observed peaks performed with 
the program UnitCell [131 . For each pressure, about 100 
peak positions were refined, covering d-spacings down to 
1.067 A. The lattice constants are reported in table H] 

We followed carefully the width of the diffraction peaks 
under pressure. Some peaks showed a marked broaden- 
ing between 40 and 60 GPa, as shown in figure [5] (top) , 
associated to a very slight anomaly in the evolution of the 
volume and lattice parameters. Alhough it is tempting to 
interpret this behaviour as the signature of a subtle phase 
change, not uncommon in perovskites, we note that this 
change is not associated to any change neither in the 
diffraction pattern (in terms of emergence/ vanishing of 
diffraction lines), nor in the Raman spectra as will be 
detailed in the next section. After a comparison with a 
second data set recorded with neon as a pressure trans- 
mitting medium (not shown here), we attributed this 
change to the onset of non-hydrostatic stress, similar to 
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Figure 2: (top) Width of Bragg peaks, averaged over selected 
reflections, (bottom) Evolution of the relative volume with 
pressure, together with the equation of state calculated by 
Ross [101 (dotted line) and from this work (solid line). The 
vertical line marks the onset of non hydrostatic stress. Only 
data below this limit were used for the determination of the 
equation of state. 



the observations reported by Downs et al. on monoclinic 
aegirine (28| . Details of this analysis and information 
that can be learned about the stress anisotropy in the 
cell will be detailed elsewhere. For the purpose of this 
work, we shall retain that (i) the Pnma structure is sta- 
ble up to 60 GPa and (ii) the pressure can be considered 
hydrostatic up to 40 GPa only for this experiment. 

The evolution of the volume with pressure at room 
temperature is reported in figure [5J In the so-called F-f 
plot [2§| (not shown) , our data between and 40 GPa can 
be represented by a straight horizontal line, indicating 
that the data can be satisfactorily fitted by a second- 
order Birch-Murnaghan equation of state (EoS). The fit 
was performed by the program EoSFit js^l- This fit yields 
a bulk mudulus of Kq = 181.0(6) GPa for a pseudo-cubic 
volume Vo = 55.830(17) A"^ and a constrained value of 
i^Q = 4. The same procedure was used to determine the 
compressibilities of the three axes. Similarly, a second- 
order Birch-Murnaghan EoS was found sufficient in all 
cases. We find Ka = 185.8(1.0) GPa, Kb = 182.6(1.1) 
GPa and = 174.5(2.3) GPa. 

In table im we compare our parameters of the EoS with 
previous results. The values reported by Ross and Angel 
|ig | were recorded on a single crystal up to 10 GPa. They 
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Table II: Parameters of the equation of state (EoS) from this 
work compared to the values determined by Ross and Angel 
[3, Xiong et al. ^ and Wu et al. The EoS is indicated 

as follows: M = Murnaghan EoS, BMi = i*''-order Birch- 
Murnaghan EoS. Values in italic indicate parameters that are 
either unrefined or implied by the model used. 



have found Vq = 55.941(4) A^, Kq = 170.9(14) GPa, and 
K' = 6.6(3) from the fit of the P-V data with a third- 
order Birch-Murnaghan EoS. Although this result is in 
very good agreement with previous determination of Kq 
by ultrasonic measurements (see Ref. [l^ and references 
therein), this equation of state extrapolated up to 40 GPa 
does not reproduce our experimental data at high pres- 
sures (flgure[2|). The fit to our data yields different values 
with a bulk modulus Kq exceeding by 10 GPa (6%) the 
value given by Ross. This difference might be explained 
by the fact that they used P-V data up to 10 GPa to con- 
strain their EoS, whereas we used values up to 40 GPa. 
A fit of our data up to 10 GPa only yields Kq = 172(3) 
GPa and K'q = 5.1(7), in agreement with the results by 
Ross. On the other hand, our values are far from Xiong 
et al. [l^l, probably because of the relative inaccuracy 
of the lattice parameters derived from their experiment 
in non-hydrostatic conditions. The same limitations ap- 
ply to the result by Wu et al. [l^, obtained from a fit 
to a Murnaghan equation of state, which is known to be 
inappropriate for large compressions. 

The values for the axial compressibilities are also re- 
ported in table HI] and compared to the values by Ross 
and Angel. All values are slightly higher than their re- 
sults, consistent with the higher bulk modulus. We find 
however that the least compressible of the three axes is 
the a axis, and not the b axis as reported in (T9| . This 
disagreement should be taken with care, given the fact 
that the degree of anisotropic compression remains very 
small. This may have an influence on the estimation of 
the pressure- induced distortion, as will be discussed later. 

For the purpose of estimating the tilt angles from 
diffraction data, it is useful to recall how the tilt sys- 
tem in Pnma perovskites can be described: rather than 
with Glazer's tilt system, two tilt angles <f> and 9 are of- 
ten used, representing tilts around the [010] and [101] 
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Figure 3: (top) Pressure dependence of spontaneous strains 
under pressure, (bottom) Evolution of the ai angles under 
pressure, as calculated from equation The dotted line is 
a linear fit to the data. 



directions respectively, relative to the pseudo-cubic axes, 
as defined by Zhao [^, [sil, [32|. Alternatively, the two 
symmetry-independent Ti-01-Ti (ai) and Ti-02-Ti (0:2) 
angles can be used. Experimentally, tilt angles are best 
determined from atomic positions obtained by full struc- 
tural refinements. Under the assumption of undistorted 
octahedra, these angles can be calculated from the lattice 
parameters only, using the following geometrical relations 
(Pnma settings): 



coscj) = cspljh 
cos 9 = c/a 

The ai angles can in turn be calculated from Q: 

cos9 = cos{(180- ai)/2} 

COS0 = COs{(180-Q!2)/2}/\/ cos 



(1) 



(2) 



With this parametrization, an increase of the distortion 
corresponds to a decrease of the angles ai. At ambi- 
ent conditions, the angles calculated from the atomic po- 
sitions found in Ref. [2^ are cj) = 8.9 ° , 9 = 11.6 ° , 
ai = 156.8 ° and a2 = 155.8 ° . 

We calculated the tilt angles from the experimental lat- 
tice constants using the geometrical equations © . Those 
relations are known to underestimate the actual tilt an- 
gles determined from atomic positions [5|, so we addi- 
tionally rescaled the values so as to match the angles at 
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room pressure. The evolution is shown in figure [3] (bot- 
tom). The angles decrease slightly with pressure, with 
an average rate of -0.045 ° /GPa. 

It is also possible to quantify the distortion in Pnma 
perovskites by spontaneous strains. As pointed out in 
Ref. [l^, [33, [s^l, spontaneous strains of interest 
for the orthorhombic structure can be calculated from 
the pseudo-cubic lattice parameters as follows: £ct = 
|c — ao|/ao measures the distortion from cubic to tetrago- 
nal and Eot = \0' — h\/aQ the distortion from orthorhombic 
to tetragonal, were oq is the cubic lattice constant ap- 
proximated here as V^/^ . These quantities are plotted in 
figure |3] (top) . Both strains increase up to 40 GPa, but 
so slightly that it almost remains constant within the ex- 
perimental uncertainty. This small variation is consistent 
with the observation by Ross and Angel [l^ who had ob- 
served no pressure-dependence of these strains between 
and 10 GPa, and confirms that these strains are not well 
adapted to the study of the very weak pressure-induced 
distortion in CTO. 



IV. RAMAN SPECTROSCOPY 

The group theoretical analysis for CTO in the Pnma 
space group yields 24 Raman active phonon modes de- 
composed as 7 + + 7 B2g + hBj,g. The Raman spec- 
trum of CTO has been studied experimentally by several 
authors in the past at ambient conditions, as a function 
of temperature (35[437| , under pressure up to 26 GPa (2l| 
and under high-pressure and temperature js^. In spite 
of these works, a comprehensive assignment of the 24 
Raman active phonons is still lacking. This is due to (i) 
band overlapping of the numerous Raman modes, (ii) the 
small distortion from the cubic perovskite which may re- 
sult in narrow splittings and low intensities for some Ra- 
man modes and (iii) an intense background scattering. 
This background can be described as two broad bands in 
the 200-550 and 600-750 cm~^ ranges, usually accounted 
for by second-order scattering processes, by analogy with 
the cubic perovskite SrTiOa in which a very similar back- 
ground is observed although first-order Raman modes are 
forbidden by symmetry. A partial assignment of the 7 Ag 
and some i?2g modes, recalled in table IIIll was proposed 
by Mc Millan et al. [s^. The frequencies and symme- 
tries of Raman active phonon modes were also predicted 
by ab-initio calculations [s^ but the results appear to be 
strongly dependent on the choices made on the structure 
relaxation (atomic positions and/or lattice parameters) 
in the calculation, and has not been used so far to clarify 
the assignment. We also note that both studies contain 
some confusion in the group-theoretical analyses due to 
the choice of the orthorhombic space group setting Pnma 
or Pbnm. A correspondence between the two settings 
and the details of the Raman tensors in both cases can 
be found for example in (40| . The Pnma setting is used 
throughout this paper. 

Selected Raman spectra recorded under pressure are 
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Figure 4: (top) Polarized spectra recorded at 54.7 GPa in par- 
allel (VV) and crossed (HV) polarization. The modes iden- 
tified as Ag (resp. Big) are marked with stars *. (bottom) 
Selected unpolarized spectra recorded under pressure. The 
spectrum at ambient conditions was recorded on a different 
crystal. 



presented in figureU) The intensity of the broad bands as- 
sociated with second-order processes gradually decreases 
with pressure. At 30 GPa, the band originally situated in 
the 200-550 cm^^ range can be considered to have com- 
pletely vanished, while some traces of the second band are 
still visible in the high-frequency part of the spectrum at 
800-1000 cm~^. This phenomenon is reminiscent of the 
evolution of the Raman spectrum of SrTiOjunder pres- 
sure and probably has the same origin (i,|4lj. A total 
of 12 sharp first-order peaks are identified, as well as a 
broader band at higher frequency that was fitted with 
two peaks, although more contributions might in fact be 
overlapping. The evolution of the Raman shifts for the 
identified bands is shown in figure [5] All modes harden 
under pressure, although with somewhat different slopes 
and curvatures. We determined the initial slopes of the 
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Figure 5: Evolution of tiie Raman shifts with pressure. The 
strongest peaks are marked with full symbols while the open 
symbols indicate weaker bands. The purple diamonds corre- 
spond to second-order features. The dashed line are linear 
fits to the data in the 0-25 GPa range. 



different modes by linear fits in the 0-25 GPa range (table 
IIII|) . The linewidths (not shown) vary little with pressure 
and do not show any change that would indicate a phase 
transition. We measured the spectrum after the exper- 
iment (after opening of the cell) and ensured that the 
typical spectrum of CTO was recovered. 

Polarized spectra were recorded at 55 GPa, the highest 
pressure investigated (figure 2]) . The orientation of the 
crystal used in the Raman experiment was previously de- 
termined from x-ray diffraction in the DAC at room con- 
ditions. In this geometry the incident laser propagation 
direction lies very close to the z-axis (less than 10 ° off) so 
that the polarization of the laser lies in the {xy) plane, 
allowing us to separate Ag and Big modes in parallel 
(VV) and crossed (HV) polarizations respectively. The 
spectra in figure S] were obtained after rotating the DAC 
so as to maximize the contrast between the two configu- 
rations, in order to be as close as possible to the crystal- 
lographic axes. Our polarized spectra therefore allow us 
to unambiguously assign the modes at 183, 227, 344 and 
489 cm~^ to Big symmetry, all of them being strong in 
the (HV) spectrum and almost absent in the (VV) spec- 
trum (all frequencies are refered to at room conditions). 



Similarly, the modes at 155, 251, 328 and 467 cm~^ can 
be assigned to Ag symmetry. Those four peaks are also 
present in the (HV) spectrum, probably as a result of po- 
larization leakage and imperfect sample orientation. The 
situation is less straightforward for the mode at 288 and 
the intense modes at 500-510 cm^^ which are observed 
in both geometries with almost the same intensity. We 
assign the mode at 288 cm~^ to Ag symmetry, although 
its rather strong intensity in the (HV) spectrum might 
indicate a mode superposition. The modes at 500-510 
cm"^ are peculiar in so far as they are hardly visible at 
room conditions but gain intensity under pressure. These 
modes correspond to the strong broad band at 493 cm~^ 
in the Big spectrum noted by Mc Millan and Ross (36| . 
who mentioned that it could be due to either first- or 
second-order scattering. We believe that a second-order 
nature of these two modes can be ruled out. First, they 
exhibit a high intensity at both high pressure and low 
temperature (37l |. Second, they are observed in both po- 
larization geometries (VV) and (VH) thus they cannot be 
pure overtones. We therefore explain them by a superpo- 
sition of modes with Ag and Big symmetry. In addition, 
three very weak bands are observed at 144, 162 and 453 
cm~^. Tlreir weak intensity does not enable their sym- 
metry to be determined from the polarized spectra. It 
is nonetheless possible to obtain some information from 
the fact that the mode at 144 cm~^ crosses the mode at 
155, which was assigned Ag symmetry. Since two modes 
crossing each other have to be of different symmetries, 
the mode at 144 cm~^ is expected to be Big. Similarly, 
the weak mode at 162 cm~^ crosses the mode at 183 
cm~^. This suggests B2g or B^g symmetries for these 
weak modes, only observed because of imperfections in 
the sample orientation and polarization conditions. 

These assignments are summarized in table lllll and 
compared to the previous assignment by Mc Millan and 
Ross [s^ . We do not confirm the assignment of the modes 
at 183 and 227 cm~^ to Ag symmetry as they proposed, 
but we note that they have not observed these two modes 
in our configuration {z{xx)z or z{yy)z). 

In a second step, it is useful to associate to each Ra- 
man mode specific atomic displacements. A significant 
amount of work has been carried out in the past for 
Pnma perovskites and vibrational patterns of the Ra- 
man active phonon modes have been proposed by Iliev et 
al. [i^, |4^ in their study of various AMnOs perovskites. 
The numbering of the Raman modes in the following is 
taken from these studies. This work requires the knowl- 
edge of the phonon eigenvectors, which can be obtained 
by either ab-initio or empirical lattice dynamics calcu- 
lations, often carried out using a shell model. As far 
as CTO is concerned, several ab-initio (39l . 143 . [45| and 
shell model calculations |2ll . |46| have been reported. Un- 
fortunately, these calculations have been performed for 
other purposes than mode assignment and show too poor 
an agreement between calculated and experimental Ra- 
man frequencies to be used for our study. Information 
can nonetheless be gained from comparison with Raman 
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data and simulations of isostructural CaBOa perovskites. 
Several such perovksites have been studied in the past for 
B = Ge, Sn, Zr, Mn [H, S^-SS- All of them have the 
Pnma orthorhombic structure and the distortion from 
the cubic perovskite in CTO, as measured for example 
by the tilt angles, is intermediate between CaGeOa and 
CaMnOa. Since the B cation sits on a center of symme- 
try of the crystal, it does not participate in any Raman 
active mode, and can be expected to have a moderate in- 
fluence on the Raman frequencies. This is supported by 
the Raman work by Tarrida et al. [4§| across the solid so- 
lution CaZrOa-CaSnOa, which shows that most observed 
Raman frequencies are only weakly dependent on the B 
cation. Conversely, their analysis of the Raman spectra 
across the CaZrOa-SrZrOa system clearly points out the 
Raman modes that strongly depends on the A cation. 
The comparison for the different CaBOs perovskites is 
shown in table IIVI The agreement between the experi- 
mental Ag frequencies is remarkable. The agreement is 
much more difficult to verify for the Big modes due to 
the incomplete mode identification. It is nonetheless pos- 
sible to proceed to an identification of the pattern of vi- 
brations from the CaMnOa [i^. The calculated Raman 
frequencies for CaMnOa, on which the mode assignment 
is based, is also recalled. For comparison, the results of 
the shell model calculations for CaMnOa by Sopracase et 
al. (131 are also indicated. Their model was successfully 
fitted against the infrared active frequencies, but shows 
marked differences to the Raman frequencies. 

Let us consider first the set of modes that involve dis- 
placements of the Ca cation in the xz plane or along the y 
direction. The three Ag modes Ag{5,6,7) can be identi- 
fied in all the compounds at very similar frequencies that 
fit well with the shell model calculations by Abrashev 
et al. on CaMnOa Among the Big modes involv- 

ing Ca-displacements, few comparisons can be made, ex- 
cept for i?2g(7) that has been observed in CaGeOa and 
CaZrOa. This attribution is very well corroborated by 
the study of (Ca,Sr)ZrOa system, which shows that these 
four modes are by far the most A-cation dependent (49j . 
In CaTiOa the corresponding mode frequencies are 251, 
288 and 328 cm^^. In addition, the associated mode at 
344 cm~^ matches very well the Big{5) mode that also 
implies Ca shifts along the y direction, which supports 
our symmetry-based assignment. 

We now turn to the highest frequency modes that in- 
volve B-0 bond stretching. Three such modes are ex- 
pected in the 500-550 cm"-'^ range. The ^4^(1) mode was 
observed in CaTiOa, CaZrOa and SrZrOa and it was 
shown [4^ that its frequency was almost insensitive to 
the A cation. This mode is closely related to the Big{2) 
mode that involves the same movements. This supports 
our assignment of a superposition of Ag + Big for the 
bands at 505 cm~^. Other modes involving bond stretch- 
ing are expected from the calculation at higher frequen- 
cies {Big{l), B2g{l), B3g{l)) but have not been observed 
conclusively in any of the compounds. 

Of particular interest are the modes corresponding to 



Wo Symmetry Slope 



(cm-i) 


This work 


Ref. [36| (cm 


"VGPa) 


144 


B29 or Bsg 




1.72 


155 


A, 


Ag 


1.28 


162 


B2g or Bsg 




1.62 


183 


Big 


Ag{ + B2gl) 


1.35 


227 


Big 


Ag{ + B2gl) 


2.12 


251 


A, 


Ag{ + B2g1) 


1.97 


288 


Ag 


Ag 


1.64 


328 


Ag 


Ag 


3.11 


344 


Big 




3.27 


453 


B2g or Bzg 




1.89 


467 


Ag 


Ag 


2.36 


489 


Big 




2.56 


500 


Ag 




4.84 


510 


Big 




4.96 




Second-order bands 




638 






5.04 


682 






4.95 


753 






5.78 



Table III: Summary of the observed Raman modes. 



rotations of the octahedra around the [010] and [101] di- 
rections respectively: Ag{2) and Ag{A). A comparison 
with the modes in CaMnOa suggests they correspond to 
the two low-frequency Ag modes. One of them in CTO 
can be identifed as the mode at 155 cm~^. An ambiguity 
remains for the second rotational mode: the compari- 
son with the other compounds suggests this could be the 
mode at 188 cm~^ identified as an Ag mode by McMil- 
lan and Ross j36|, even though we could not confirm this 
assignment. The two low-frequency i?ig(3) and i?ig(4) 
modes also involve rotations of the octahedra and our 
values for CTO compare well with the frequency mea- 
sured in CaGeOa and CaMirOa- 

V. DISCUSSION 

A. Stability of the Pnma structure 

Our Raman and X-ray scattering experiments show 
that the Pnma orthorhombic structure of CTO remains 
stable up to 60 GPa. This result is in contradiction with 
the report by Xiong et al. of a phase transition to a 
hexagonal structure at 10 GPa, but confirms and extends 
the study by Gillet et al. [2l| whose Raman scattering 
experiment had not revealed any phase transition up to 
26 GPa. The change interpreted by Xiong et al. [22| as 
a phase transition was most probably due to the solid- 
ification of the ethanol-methanol pressure transmitting 
medium around 10 GPa, all the more that they could 
not reproduce this transition in a subsequent experiment 
carried out without any transmitting medium. Moreover, 
the absence of a phase change does not confirm the ab- 
initio calculation by Wu et al. who had predicted 
a transition to the post-perovskite Cmcm phase at 30 



8 





CaGeOa CaTiOs 




CaMnOa 


CaSnOa 


CaZrOa 




Mode 


Ref. [47,1 This work 


Ref. 


[43J 


Ref. [50| 


Ref. [49| 


Ref. [51| 




symmetry 


Exp. Exp. 


Exp. 


Calc. 


Calc. 


Exp. 


Exp. 


Mode description 




151 155 


160 


154 


146 


145 


145 


In-phase j/-rotation 




171 


184 


200 


229 


- 


- 


In-phase i-rotation 




265 251 


243 


242 


323 


265 


263 


Ca and 01, a::-shifts 


M^) 


284 288 


278 


299 


336 


277 


287 


Ca, 2-shifts 




328 328 


322 


345 


406 


355 


358 


Ca and 01, a;-shifts 




467 


487 


467 


529 


442 


439 


Out-of-phase B-0 bending 


^9(1) 


500 


- 


555 


657 


- 


543 


In-phase B-O stretching 


Sl9(3) 


184 183 


179 


178 


207 


- 


- 


Out-of-phase {/-rotation 


Sl9(4) 


237 227 


- 


281 


265 


- 


227 


In-phase x-rotation 




358 344 


- 


354 


450 


- 


305 


Ca and Ol, y-shift 




510 


- 


536 


656 


- 


439 


Out-of-phase B-O stretching 


Bi,m 


489 


- 


743 


700 


- 


547 


Out-of-phase B-O stretching 


^29(4) 


162 


- 


148 


229 


- 


190 


Out-of-phase ^-rotation 


^29(7) 


248 


- 


232 


249 


- 


212 


Ca and 01, z-shifts 


B2A5) 


- 


- 


292 


306 


- 


234 


Ca and 01, z-shifts 


B2g{6) 


- 


- 


366 


420 


- 


418 


Ca, a::-shifts 


^2,(3) 


- 


- 


453 


503 


- 


469 


Out-of-phase B-O bending 


B2A2) 


- 


- 


485 


541 


- 


543 


In-phase B-O bending 


^2,(1) 


- 


- 


749 


713 


- 


- 


In-phase B-O stretching 


^39(5) 


- 


- 


290 


234 


- 


- 


Ca, y-shifts 


^39(4) 




320 


304 


365 






In-phase ^-rotation 


^39(3) 


380 




459 


503 






Out-of-phase B-O bending 


Bs^m 


496 


564 


541 


655 






Out-of-phase B-O stretching 


S3s(l) 






754 


719 






Out-of-phasc B-O stretching 



Table IV: Comparison of the observed Raman modes for different CaBOa perovskites, ordered roughly from the less distorted 
to the most distorted structure. The numbering and description of the Raman modes follows Ref. |43l |. The assignment of 
many Big Raman modes is unclear due to incomplete information. See the text and references for details. 



GPa [2J|. This disagreement could be partially due to 
the ab-initio calculations being performed at K while 
our measurements are performed at room temperature 
and do no exclude the possibility of a phase transition 
at low temperatures. Alternatively, it could be that high 
temperatures are necessary to break the bonds involved 
in the reconstructive perovskite to post-perovskite tran- 
sition [2!. 



B. Evolution of the distortion under pressure 



We now want to discuss the evolution of the distor- 
tion with pressure, especially in the high-pressure regime. 
First, we shall recall briefly the main features of the va- 
lence bond sum model by Zhao [l6| . For a given site 



A or B, the valence bond sum Vi is 



V, 



N 

E 



exp 



R(] — Rj 



(3) 



where N is the coordination number, Rij is the bond 
length, i?o a constant specific to an anion-cation pair 
and b = 0.37 A a universal constant. The valence bond 
sum equals the formal valence of the cation at ambient 
conditions. It is then assumed that the pressure-induced 
changes in the valence bond sum for sites A and B are 
equal ("bond valence sum matching principle"). It can 
then be shown that the compressibility ratio Pb/Pa is 
estimated by Ma/Mb, where Mi is the site parameter 
defined from the average bond lengths Ri by 



M, = 



R,N 



■ exp 



Rq — Ri 



(4) 



b \ b 

In addition, it has been shown later [l^ that the same 
model could be used to calculate the evolution of the 
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tilt angle under pressure from the knowledge of the de- 
tailed structure at ambient conditions and the pressure- 
dependence of the lattice parameters. For A^^ B'^'^O^ 
perovskites such as CaTiOa, this model yields Pb/Pa < 
1, i.e. the distortion should increase under pressure, 
and the calculated increasing rate of the tilt angle is - 
0.131 ° /GPa, in good agreement with experimental data 
collected up to 10 GPa [15„ .16]. In the following, we will 
analyze the evolution of the distortion from our data, and 
check this model against the evolution up to the highest 
pressures measured. 

Within the model by Zhao, the angles ai are deter- 
mined from the linear compressibilities of the lattice pa- 
rameters Pa, Pb, Pc and the Ti-0 bonds, using the rela- 
tion 



a, = 2sin-^[exp(A/3,P)sin(ao,i/2)] 



(5) 



where A/3i = Pb~oi - Pb and AP2 = Pb-02 - (a'^Pa + 
c^Pc)/{a? + c^). The average bond compressibility Pb is 
a good estimate for Pb-Oi- Pb in turn is calculated from 
the variation of the valence bond sum Vb under pressure: 



Pi 



b dlnVB 
dP 



(6) 



Finally, Vb can be estimated at a given pressure in the 
following way: we calculate for each pressure the average 
of the bond valence sums for sites A and B using the 
"fixed-coordinates" model, whereby the atomic coordi- 
nates in the cell are kept equal to their values at ambient 
conditions, and only the lattice parameters are set to 
their experimental values. It was shown in [l^ that this 
yields a good estimate for the pressure-induced change in 
the bond valence sum AVb- 

For the calculation of the tilt angles, we take the initial 
structure from Sasaki [2^ and the values for Rq (1.815 A 
for Ti-0 and 1.967 A for Ca-0) are taken from Ref. 
For the calculation of the compressibilities, we proceeded 
in two steps. In a first step, for comparison purposes 
with the results by Zhao uM, '^^ estimated the linear 
compressibility of axes as well as d In / dP by a linear 
fit to our data between and 10 GPa. With pressure- 
independent compressibilities calculated in this way, we 
find an average rate of d{ai)/dP = —0.134 °/GPa. in 
excellent agreement with the -0.131 ° /GPa by Zhao [l6| . 
showing that the differences found in the previous section 
regarding the axis compressibilities are small enough to 
be neglected for the calculation of the tilt angles. The 
result of this first calculation is shown in figure [5] (solid 
line). The comparison with the result found from the 
direct calculation of the tilt angles from the lattice con- 
stants alone (dotted line in figure 6) shows that the geo- 
metrical relations ^ do not yield a correct angle change 
rate (-0.045 vs. -0.134 ° /GPa). This resuh most proba- 
bly reflects that the hypothesis of undistorted octahedra 
is invalid at high pressures. 

In a second step, we tentatively performed the same 
calculation over the full pressure range with pressure- 
dependent axial compressibilities. Of course, we expect 




20 

Pressure (GPa) 



40 



Figure 6: Evolution of the average Ti-O-Ti angle under pres- 
sure calculated in different ways: calculated from the lattice 
constants using the geometrical relations ^ (dotted line), 
calculated from the bond valence sum model with pressure- 
independent compressibilities (solid line) , calculated from the 
valence bond sum model with pressure-dependent compress- 
ibilities (dashed line), calculated from the frequency of the 
Ag{2) Raman mode (symbols). 



the model to break down at some point, since the pa- 
rameters b and Rq used to calculate the bond valence 
sum are determined by the requirement that the valence 
bond sum has to be equal to the formal charge of the 
cation at ambient pressure. The result is shown in figure 
E] (dashed line) . The angles follow closely the results of 
the first calculation close to ambient pressure but show 
a much more pronounced curvature and even start to in- 
crease again around 35 GPa. This is not unrealistic, since 
a change in the compression mechanism could in princi- 
ple result in such an evolution, as is presumably the case 
in ScZrOafl^l. In CTO, this behavior remains to be ver- 
ified experimentally, which would require full structural 
refinements to calculate tilt angles as accurately as pos- 
sible. Unfortunately, full structural refinements on single 
crystals in this pressure range are technically challenging, 
and X-ray diffraction is not a very sensitive tool for the 
determination of oxygen positions. 

In this context, Raman spectroscopy may be an inter- 
esting alternative to follow the two compression mecha- 
nisms: bond compression and octahedra tilting. In prin- 
ciple, it is possible to estimate directly the compressibil- 
ity of a polyhedron from high-pressure Raman measure- 
ments provided that a phonon mode can be found that 
involves in-phase expansion of all its bonds j52| . The 
relevant quantity is the mode compressibility that is pro- 
portional to (dlnw/dP). By comparing the relative in- 
crease of the chosen phonon frequencies, it is possible to 
estimate a compressibility ratio and decide whether the 
distortion increases or decreases under pressure. How- 
ever, in a crystal with distorted polyhedra where most 
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phonon modes involve mixed displacement of atoms, it 
is not straightforward that a suitable phonon mode ex- 
ists that faithfully reflects the compression of the bonds 
or the polyhedra. As far as tilt angles are concerned, 
Iliev et al. have shown in their studies of a series of 
rare-earth manganates that the two rotational Ag{2) and 
Ag (4) mode frequencies scale with the tilt angles with the 
empirical rate 23.5 cm~-^ / ° |53|. This rule has also been 
found approximately valid in rare-earth scandates (s^ . 
Unfortunately, A^+Ti^+Os perosvkites with the Pnma 
structure are rare and the same systematic work can 
hardly be done. We may however expect a similar rule to 
hold true for CTO, i.e. the frequencies of the rotational 
modes to scale with the tilt angles. 

From our mode assignment, we oberve that the most 
pressure-dependent modes are the modes involving Ti- 
O bond stretching and Ca-shifts. On the other hand, 
the rotational modes at low frequencies show a much 
more moderate frequency increase under pressure. This 
is qualitatively consistent with the general picture of a 
volume reduction obtained by an overall compression of 
the bonds with only a slight variation of the tilt an- 
gles. The bands with the highest mode compressibility 
are the bands at 500-510 cm~^ which we have assigned 
to Ti-0 stretching. According to the vibrational pat- 
tern by Iliev. the mode best suited to the examination of 
the compression of TiOg octahadra would be the 53^(1) 
at high frequencies (octahedra "breathing" mode). This 
mode however, was not identified in CaTiOa nor in any 
of the CaiJOa compounds mentioned above. Moreover, 
the mode compressibility calculated from the evolution 
of bands at 500-510 cm~^ is the strongest of all observed 
modes, and we cannot expect to extract compressibility 
ratios in agreement with the value now well established 
from crystallography. This reflects the need for a more 
precise identification of the relevant breathing modes. 

Finally, we consider the evolution of the rotational 
modes and assume that the Ag{2) and Ag{A) scale with 
the tilt angles (j) and respectively. In CTO, we have 
(j) = 8.9 ° and 6 = 11.6 ° . The Ag{2) rotation mode was 
identified as the mode at 155 cm~^ and we additionally 
assume the Ag{4) is indeed the mode at 181, although it 
is emphasized that some ambiguity remains. From the 
evolution of their frequencies under pressure (1.28 and 
1.35 cm~^/GPa respectively), we calculate an average 
angle variation of -0.160 ° /GPa, which compares reason- 
ably well to the -0.134 ° /GPa calculated by the valence 



bond sum model. It is therefore not unreasonable to con- 
sider that these phonon frequencies provide a good ap- 
proximation of the average tilt angle. The average angle 
calculated in this way (and shifted to match the value at 
zero pressure) is plotted in figure [B] together with the re- 
sults of the different calculations. The evolution suggests 
that the increase of the tilt angles slows down as pressure 
increases, which is consistent with the prediction of the 
pressure-dependent valence bond sum model, although 
the pressure ranges involved are very different. Further 
work and more conclusive mode assignments are needed 
for a validation of this approach, but we believe that Ra- 
man spectroscopy can be a valuable tool to estimate tilt 
angles in the high-pressure region. 



VI. CONCLUSION 

We have performed X-ray diffraction and Raman spec- 
troscopy on single crystal CaTiOa in a diamond-anvil cell 
up to 60 GPa and 55 GPa respectively. The orthorhombic 
Pnma structure remains stable over the whole pressure 
range investigated, which contradicts the ab-initio pre- 
diction of a transition to a post-perovskite phase at 30 
GPa. We calculated the bulk modulus of CTO Kq by a 
fit to the compression data in the 0-40 GPa range and 
found a value of 181.0(6), in good agreement with previ- 
ous results. In addition, polarized Raman spectra were 
recorded at high pressure and have enabled us to clarify 
the Raman mode assignment of CTO, and identify reli- 
ably Raman modes associated to Ti-0 bond stretching, 
Ca shifts and octahedra rotations. The analysis of the ro- 
tational Raman modes as well as calculations using the 
valence bond sum model both suggest that the weak dis- 
tortion of the Pnma structure increases under pressure 
but with a decreasing rate. 
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